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Abstract 
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to prove replica symmetry, we show that the fluctuations of the overlaps 
and of the free energy are Gaussian, on the scale l/y/N, for N large. 
The method we employ is based on the idea, we recently developed, of 
introducing quadratic coupling between two replicas. The proof makes use 
of the cavity equations and of concentration of measure inequalities for the 
free energy. 
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1 Introduction 



We consider the mean field spin glass model introduced by Sherrington and Kirk- 
patrick in [1], [2], in the regime of high temperature or, equivalently, of large 
magnetic field. On physical grounds, it is known that in this region the replica 
symmetric solution holds, as shown for example in [3], and references quoted 
there. However, due to the very large fluctuations present in the model, it is dif- 
ficult to give a mathematically rigorous description of this region. Rigorous works 
on this subject include [4], [5], [6], [7], [8]. For other rigorous results concerning 
the model, we refer to [9], [10], [11], [12], [13]. 

The method developed in [8] by Michel Talagrand is particularly interesting. 
The starting point is the very deep physical idea that the phenomenon of replica 
symmetry breaking can be understood by studying the properties of the model 
under the application of auxiliary interactions, which explicitly break replica 
symmetry. In [8], this idea is employed to prove that the replica symmetric 
solution holds in a region, which (probably) coincides with that found in the 
theoretical physics literature [3], i.e., up to the Almeida-Thouless critical line. 

Recently [14] we proposed a different strategy, which consists in coupling two 
replicas of the system by means of a term proportional to the square of the de- 
viation of the overlap from its replica symmetric value. In this way, we proved 
that replica symmetry holds in a region above the Almeida-Thouless line. In the 
same region, we obtained a control of the two-replica system, provided that the 
coupling parameter is small enough, and we showed that the fluctuations of the 
overlap are at most of order 1/y/N. In the present paper we prove that, in the 
same region of parameters, the fluctuations of overlaps and free energy, when 
suitably rescaled, have a Gaussian distribution when N — > oo. The main ingredi- 
ents of the proof are the control of the thermodynamic limit obtained in [14] and 
concentration of measure techniques inspired by Talagrand's works. Then, by 
means of the cavity method, one can write self-consistent linear equations for the 
characteristic functions of the fluctuation variables, which can be easily solved. 

Previous results concerning limit theorems for fluctuations in the high tem- 
perature region of mean field spin glass models include [4], [6], [15], [9], [16]. 

This work is organized as follows: In Section 2, we recall the main definitions of 
the model and introduce the overlap distribution structure. In Section 3, we state 
the main results. Two useful tools, i.e., exponential inequalities and the cavity 
method, are briefly outlined in Sections 4 and 5. In Sections 6 and 7, we prove 
the central limit theorem for overlap and free energy fluctuations, respectively. 
Finally, Section 8 is dedicated to a short outlook about open problems and further 
developments. 
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2 The model 



The generic configuration of the Sherrington-Kirkpatrick (SK) model is deter- 
mined by the N Ising variables Oi = ±1, % = 1,2, ... ,N, and the Hamiltonian 
is 

H N (a,h;J) = — — ^ JijCTjCTj - h^aj, (1) 

where the sum runs over all the N(N — l)/2 distinct couples of sites. The 

Jjj's (quenched noise) are independent centered unit Gaussian variables Af(0, 1). 
The first term in (1) is a long range random two body interaction, while the 
second represents the interaction with a fixed external magnetic field h. For a 
given temperature 1/(5 we can introduce the disorder dependent partition function 

Z N {(5, h;J) = J2 exp(-(3H N (a, h; J)) (2) 

and the auxiliary function 

a N ((5,h) = -^E\nZ N ((5,h;J), (3) 

where E denotes the average with respect to the external noise J. Note that 
aiN((3,h) is the quenched average of the free energy per spin, apart from the 
multiplicative factor —1/(5. 

For later convenience, it is useful to generalize the model (1) by introducing a 
quenched random external magnetic field, which at every site is an independent 
Gaussian variable of strength x > 0. In other words, the Boltzmann factor of the 
system becomes 




where Jj are i.i.d. jV(0, 1) random variables, independent of the J^s. We let 
t = (3 2 in the two-body term. In the following, we always regard the system as 
depending on the parameters t, x, (5h. In analogy with Eqs. (2), (3), we define the 
disorder dependent partition function Z^{t,x, h; J) and the auxiliary function 

a N (t, x, h) = —E In Z N (t, x, h; J). 

Here, E denotes averaging with respect to and J«. For simplicity of notations, 
here and in the following we write the argument h instead of (3h. 

Let us consider a countably infinite number of independent copies (replicas) 
of the system, whose spin variables erf are distributed, for fixed J, according to 
the product state 

= ^j^j • • • , 
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where Qj = Q^^' x ' h denotes the Gibbs state associated to the Boltzmann factor 
(4). Each replica is subject to the same quenched noise. The "real replica" 
approach has already been exploited in a number of papers [17], [18], [19], [20]. 
The overlap between two replicas a, b is defined as 

i 

with the obvious bounds 

-1 < q ab < I- 

For a generic smooth function F of the overlaps, we define the (.) average as 

(F(g 12 ,g 1 3,...)) = ^ J (F(g 12 ,g 13 ,...)). 

Note that the average over disorder introduces correlations between different 
groups of replicas, which would be independent under the Boltzmann average 
Qj. For example, 

^j(<?12<?34) = ^j(<?12)^j(<?34) 

but 

(<?12<?34) + (?12>(?34>- 



3 The high temperature region and the main 
results 

In this Section, we recall the results of [14] and state limit theorems for fluctua- 
tions, in the region where we prove that replica symmetry holds, i.e., 

lim a^(t, x, h) — a(t, x, h). 

N-^oo 

a(t,x,h) is the replica-symmetric free energy [1], [2] 

a(t,x, h) — In 2 + J lncosh(/3/i + zy/tq + x) dfi(z) + ^(1 — q) 2 , 

q is the Sherrington-Kirkpatrick order parameter, defined as the unique [17] so- 
lution of 

q — q(t,x, h) = J tanh 2 (/3/i + zyjtq + x) dfx(z) 

and d/i(z) is the centered unit Gaussian measure. 

In [14] we proved the following: Consider the auxiliary function a N , dependent 
on the parameter A > 

a N (t,x, h; A) = a N (t,x,h) + lnfy^ ( e "£(«»-5) a ) 
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and the trajectory in the (t,x) plane 

r^(t',x t ') = (t',x + q(t-t')) = (t',x -qt'), 0<t'<t (5) 

where x = x + qt and q = q(t, x, h) = q(t', x t i, h). Notice that equals ajy for 
A = 0. Given x , h there exists a value t c (xo, h), such that 

k 

\a(t',x t >,h) - a N (t',x t >,h; A)| < — (6) 
for some constant k, uniformly in the triangular region 

0<t' + \<t< t c (x , h). (7) 
In the same region, the overlaps self-average around the value q: 

((ft* " Q?) < (8) 
The critical value t c (xo, h) is determined in the following way [14]: Let 

A(x , h, Ao) = ^ max ( / ln(cosh p + tanh 2 (/3/i + z^fx^) sinh p)dfi(z) — pq — 

2 peR \J ' 2A 

where A > 0. Then, we define t c (x , h) such that, for any A < t c (x , h), one has 

A(x ,h,X ) = 0. 

In the case of vanishing external field x—h—0, then also xo=q=0 and t c = 1, the 
correct critical value. As discussed in [14], the region defined by (7) falls short of 
the Almeida-Thouless line, which is the expected critical line. 

In this paper, we investigate more precisely the behavior of fluctuations of 
physical quantities around the replica symmetric value. First of all, we give a 
central limit-type theorem for the rescaled overlaps 

& = VN(qab-q), 

showing that they behave as centered Gaussian variables characterized by a non- 
diagonal correlation matrix. Notice that, thanks to (8), one has the following 
bound for the second moment of the rescaled overlap fluctuations: 

<(0 2 > < k. (9) 

Theorem 1. If t < t c (x ,h), the rescaled overlaps £^ tend in distribution, for 
N — > oo, to jointly Gaussian variables £, a b, with covariances 

(£ b ) = A(t,x,h) 
(UU) = B(t,x, h) 

{babied) = C(t,X, h), 
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where b ^ c, c ^ a,b and d ^ a,b. A, B and C are explicitly given by 

A(t,x,h) = (1 + 2R + AR 2 )Y + c R 2 (10) 
B(t,x,h) = (1+AR)RY + c R 2 (11) 
C(t,x,h) = 4R 2 Y + c R 2 , (12) 

where 

1 

Y(t,x, h) 
R(t, x, h) 



Y^-t 
d 



Yq- 1 + 2d - t 

and Y (xo, h), co(xo, h) and do(xo, h) are chosen in such a way that A,B,C satisfy 
the initial conditions 

A(0,x o ,h) = l-f 
B(0,x o ,h) = q-q 2 

C(0, x ,h) — J tanh 4 (zy^Xo + Ph) dfi(z) — q 2 . 
In particular, one has 

Y = J cosh -4 (zy/Ff+x + f3h) dfx(z) (13) 

Recently, an analogous result was proved independently by Talagrand [16], 
who computed the N — > oo limit for all moments of the £ variables. 

The expressions for A, B and C were first given by Guerra in [17]. For h = 
x — 0, the limit Gaussian variables are not correlated and have variance 1/(1 — t), 
which is a well known result [4] , [6] . 

Let us consider now free energy fluctuations. Aizenman, Lebowitz and Ruelle 
[4] proved that in the case of zero external field and t < 1, the variable 

In Zjy — InEZ^ 

tends to a shifted Gaussian random variable whose variance diverges at t = 1. In 
the general case the situation is quite different and the following theorem holds: 

Theorem 2. Let 

f N (t, x, h; J) = VN ^(M,fr;J) _ _^ ^ ^ 

Ift < t c (xo, h) then 

f N (t,x, h- J) J\f(0,a 2 (t,x, h)), 
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where 

ft 

a 2 {t, x, h) = Var (in cosh^V* q + x + (5h)) - 
ii/ere, Var(.) denotes the variance of a random variable and z = Af(0, 1). 

Notice that fluctuations of the extensive free energy InZjv are of order 1 at 
zero external field and of order \/~N otherwise. 



4 Exponential suppression of overlap fluctua- 
tions 

General arguments based on concentration of measure [9], [21], [22] show that 
the fluctuations of the free energy l/NlnZ^ around its mean value a at are ex- 
ponentially suppressed as iV grows. Indeed, one has the following [9] 



Theorem 3. For any u > 0, 

log Z N (t, x, h; J) - a N (t, x, h) 
where 

K= 1 



> u < exp(-NKu 2 ), (14) 



t + 2x 

This, in connection with the results of [14], allows to obtain a strong control 
on the fluctuations of the overlaps (we learned this nice argument in [8]): First 
of all, the same argument leading to Theorem 3 shows that 

P ( ^ In ft W» (e^fo"-*) 2 ) -a N (t,x,h,X) +a N (t,x,h) > 2u^j < exp(-NKu 2 ). 

Therefore, thanks to Eq. (6), with probability at least 1 — exp(— NKu 2 ) one has 

Q t ^ h (ei^--*) 2 ) < e' Nu+2C 
for A < A < t c (x , h) — t. Then, by Tchebyshev's inequality 

and, choosing u = Xv 2 /16, one has 

Vt, x ,h {X{\ qi2 -g\>v}) < e~ N ^ +2C , 
The estimate we are looking for easily follows: 

EOt** (X{|, 12 -,1>,}) < e- N ^ +2C + e~ NK ^. (15) 
Of course, this is much more than just self-averaging of the overlaps. 
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5 The cavity method 

The cavity method allows to express thermal averages of quantities defined on 
the iV-spin system as functions of averages on the system with N — 1 spins, at 
a slightly different temperature. This method has been widely applied both in 
the theoretical physics literature [3] and in the mathematical physics one (see, 
for instance, [5], [16], [10], [18]). 

Introduce the following definitions: 

t' = t(l- AT- 1 ) 

a a = (r) a ,e a ), rf e {-1, l}"" 1 , e a = a% = ±1 
J — Jn 

g i = J Ni , i = l,...,JV-l 

m = n&_i(0 

The cavity equations consist in the identity 

n N (f(« ,■■■,* )) = ' (16) 

where Av denotes the average over the spin variables e a and 

k 

¥ k) = exp£V (ytjNgr 1 a + ^J + fih^J . (17) 



a=l 



grj a denotes the scalar product Yl!i=i 9iVi- 

6 Limit theorem for overlap fluctuations 

To prove Theorem 1, it suffices [23] to show that for any integer s, the charac- 
teristic function 



\ <a,b) I 



(a,b) 

converges for iV — > oo to 

cf>\u) =exp|-^(Lu,u)|, (18) 

where (. , .) denotes scalar product and L is the s(s — 1)/2 x s(s — 1)/2 dimensional 
matrix of elements 

L(ab),(ab) = A(t, x, h) 

L(ab),(ac) = B(t,X, h) 
L(ab),(cd) = C(t,X, h). 



The idea of the proof is to obtain a set of closed linear differential equations for 
0tv(w), which determine uniquely the solution as (18), for N — > oo. Some of the 
calculations involved in the proof are quite long, although straightforward, and 
are therefore just sketched. 

First of all, we explain how the cavity equations (16), (17) can be simplified in 
the region where (15) holds. Following [10], we introduce some notations, letting 
Sl(.) = n% x ' h (.) and ft'(.) = ffcY). Moreover, we define 



b = n'(rj) E 
r] a = v a -b 



oAT-l 



X = y/t/Ngb + ^/xJ + ph 

k 

vJ/^expPTjV) 

a=l 

f(a\...,a k ) = f( V \e\..., V k ,e k ). 
Then, the following holds [10]: 
Theorem 4. 



EQ(f(a\...,a k )) = E-J^tf (Av f ■*<*>) (19) 



+tE l—tftAvf*W y eV^l (20) 

cosh fe X I N I 



l<a<c<k 



+t E \ — Q' | Avf ¥ k) y 

cosh fe X I ^ 



^Tj « 21 » 

,^c<fc / 

-ktE t -^a(Avf^ye^)+S (22) 
cosh fe X I N J K J 



and the "error term" S can be estimated as 

\s\<Mt,xMEn>(Av\f\(f;(^) 2 + y f^f) 

\ \a=l ^ ' l<a<c<k+2 ^ ' 

where w is a smooth function of its arguments, independent of N. 

Note that, with respect to Theorem 3.2 in [10], the last sum in the r.h.s. is 
performed on a < c instead of a < c. However, the proof of Theorem 4 proceeds 
exactly as in [10]. 

Theorem 4 is a sort of Taylor expansion of the cavity equations around r] a = b. 
This turns out to be particularly useful in the region where Eq. (15) holds, since 
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in this case f) a is small with large probability, and S vanishes for N — ■> oo, as we 
explain in the following. 

In order to prove Theorem 1, we first exploit symmetry between sites to write 



N 



U 



iVn ((<j r N °N - q) e lue 



(23) 

= i e^ N ) = iVN((a%ar - ?) (24) 

=i(&i,s + 2e iu * N ) =iVN((o$ 1 v«*-q)e i «e N ), (25) 



and then employ the cavity equations to express these quantities as functions of 
0, <f, ip themselves. For instance, apply Theorem 4 to the r.h.s. of Eq. (23) and 
consider the term arising from (19). After averaging on the dichotomic variables 
e, one is left with 

% (VN -iJ2 u *b V) E { ( tanh2 X-q)Q! exp (i v! } + (26) 

(a,b) 

- u rr ,E { (1 - q tanh 2 X)Q' exp (i u £ N ~ l ) } 

-(l-q) (u ar + u ar >)E {ttrnh 2 X VI' exp (iu'^- 1 )} 

u cd E {tanh 2 X (tanh 2 X - q) Q' exp w'^" 1 ) } + o(l), 

(c,d) c,d^r,r' 



where u' = u^J\ — 1/N. The term o(l) arises when exp(iu/\^N) is expanded 
around u=0 and the terms of order u 2 or higher are neglected. Indeed, one has 



vW(o>^ - q) (e**™ vf - 1 - i £ 

\ V (o,6) 



/ .a h 



- q) 



< 2VN 



jz (a:b) w**-*) y^(*N< - q) 

(a,b) V N 



= 0(iV- 1/2 ). 



Now, rewrite E as EE g , where E g denotes the average only with respect to the 
random variables J and g^ % — 1, . . . , N — 1, and notice that J, g,i do not appear 
in the thermal average Q'. Computation of E g (. . .) would be simpler if, instead 
of X, there were 



X = y/i/Ng b + v^J + [3h, 



where 



Of course, one has 



b = 



'Nq. 



X = z^tq + x + (3h, 
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where z is a standard unit Gaussian variable and equality holds in distribution 
so that, for instance, 

E g tanh 2 X = q. 

The idea is, therefore, to expand around X — X. As a preliminary fact, notice 
that the second moment of the random variable (b — b) is bounded uniformly in 
N. Indeed, 

E\\b -b\\ 2 = E(\\b\\ - < ^E{\\b\\ 2 - Nqf (27) 

= ^EQ'ia-^u 1 ) + 0(1/N) = 0(1), (28) 

thanks to Eq. (9). As an example, let us examine in detail the first term in (26), 
that is, 



NE E g (tanh 2 X) ft' exp (i u'^' 1 ) - % qVNE ft' exp (i u'^' 1 ) . (29) 

By a simple second order Taylor expansion and an integration by parts on the 
Gaussian noise g, one finds 

Eg tanh 2 X = E g tanh 2 X + ^(b - b)bE g <9 2 tanh 2 x\ x=x (30) 
+^E g dlt^x\ x=x+e{x _ X) (g{b-b)) 2 (31) 
= Q + ^(6 " b)(b + b)E g dl tanh 2 X (32) 



t 

4 



2N 
t 2 



\b - b\ \ 2 E g {dl tanh 2 x\ x=x+e{x _ x) - 81 tanh 2 x) (33) 



+ ^2 E 9 ^ 4 tanh 2 a ;| :c=w _ ;?) [(6-6)(6 + ^(6-6))] 2 (34) 

where < 9 < 1. Let analyse each term separately. Recalling the definitions of 
b and b, the second term in (32) equals 



2N 

t 



ft'(?f + V +2 - Nq) J d/i(z)dl tanh 2 ((3h + z^tq + x) (35) 
ft'(e + i: s+2 )(3r + 2q-2) + 0(1/N), (36) 



N 

where Yq was defined in (13). Another application of Taylor's expansion and 
integration by parts, together with Cauchy-Schwarz inequality and the fact that 
the derivatives of the function tanh 2 (x) are bounded, shows that the terms (33) 
and (34) can be bounded by 
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Therefore, using the estimate (28), the expression (29) reduces to 

it(3Y + 2q-2)EQ' [£ s+1>s+2 exp (iu'^" 1 )] + O^ 1 / 2 ), 

and 

iy/NE { (tanh 2 X - q)Q! exp (i u'^' 1 ) } = t(2q - 2 + 3Y )ip' + o(l), 

where 

The other terms in (26) are much simpler than (29), and can be dealt with in the 
same way. Finally, the whole expression (26) can be rewritten as 

t(2q-2 + 3Y )iP' -u rr ,(l - f) <j)' (37) 

-(? - f) Yl + - ( y ° - 1 1 - ^) 2 ) Mcrf + °w- 

a^r,r' (c,d) c,d^=r,r' 

The steps leading to expression (37) can be repeated with minor changes for the 
remaining terms (20) to (22). These terms, although they look more complicated 
than (19) at first sight, are actually simpler to treat, since a first (instead of 
second) order Taylor expansion around X = X is sufficient. This is due to the 
presence of terms like r] a r]b/N or rj a b/N, which are with large probability small, 
thanks to (9). Also in this case, one finds that terms (20) to (22) give quantities 
linear in 0' , 8(f)' , tp' , ip 1 , apart from terms of order o(l). As for the "error term" 
S which appears in Theorem 4, one can easily check that it vanishes in the 
thermodynamic limit. This is a consequence of the exponential decay of overlap 
fluctuations, as expressed by (15). 

Next, we show that terms like <p t N _ 1 (u') or ip^^^u') can be substituted by the 
same functions calculated at N, t, u, apart from negligible error terms. Indeed, 
for instance, 

^ N {u) = (exp^u'^-'+tuialtal-q)/^^ 
= (ex pi u'e- 1 ) t (l + o(l))=<P> + o(l). 

In the last step, we used Theorem 4 in order to substitute t with t'. Therefore, 
Eq. (23) reduces to a linear relation between 0, (p and ip, apart from a remainder 
which becomes irrelevant in the thermodynamic limit. In the same way, one sees 
that also Eqs. (24), (25) yield linear equations for <f>,(p,tp. Putting everything 
together, in the thermodynamic limit one has a set of coupled linear differential 
equations of the form 

$\u)=(i) t {u)v{u) + tM&{u) (38) 
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where <!>*(«) is the vector 

v(u) is a vector whose components are homogeneous linear functions of the vari- 
ables u, while M is a real square matrix with elements depending on q, Y Q alone. 
We do not report here the explicit expressions of v(u) and M, which are quite 
complicated. However, it is instructive to check that, for instance, the term (37) 
is in agreement with this structure. In fact, the coefficient of <f>' is a homogeneous 
linear function of the u variables, while the coefficient of ip' is linear in t and 
depends only on Y and q. As will be clear in the following, only the structure 
(38), and not the specific form of v and M, are needed to conclude the proof of 
the theorem. 

Assume at first that the matrix (1 — t M) is invertible, which in principle can 
fail only for a finite number of values of t, since M is finite dimensional. In this 
case, Eq. (38) can be reduced to a first order differential system in normal form: 

* ( (u) = ^(u)(l-tM)- 1 v(u), (39) 

which can be easily integrated. The most general solution for 0'(w), compatible 
with the initial condition 

0'(O) = 1, 

is of the form 

0'(u)=exp j-^(£u,u) + (p,u)}, (40) 

where p is some s(s — l)/2 dimensional w-independent vector, and K is a s(s — 
l)/2 x s(s — l)/2 real symmetric positive definite matrix. The symmetry and non 
negativity of K derive from the obvious property of symmetry among replicas, 
and from the bound 

< i, 

which holds for any characteristic function. The quadratic dependence on u of 
the exponent of 0*(w) stems from the linear dependence of the components of 
v(u). Clearly, Eq. (40) means that the random variables converge to some 

Gaussian process {£ a b}- Moreover, it turns out that the identification 

p = 

and 

k = i 

are straightforward. Indeed, it was shown by Guerra in [17] that, if the limit 
process is Gaussian, then it is centered and its covariance function is exactly L. 
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In order to conclude the proof, it remains to show convergence of the charac- 
teristic function for those possible values t where (1 — t M) is singular. For any 
5 > one can write 



ct> i N {u) = <p i N \u) + 5d t <p t N 



t=t-8 N S 



where < 9 N < I. After a straightforward computation one finds that 



\ \(a,6) a=l 



s(s + 1) , N , 2 
9 Us+I,s+2J 



By exploiting the uniform bound (9) and the arbitrariness of 5, one finds therefore 
that the theorem holds also for t — t . □ 



7 Fluctuations of the free energy 

In order to prove Theorem 2, we show that the characteristic function of fa 
converges to that of A/"(0, a 2 (t, x,h)), i.e., 



lim E e iu fN(.t,x,h) _ e -\a 2 {t,x,h) 
N^oo 



Define 



a{t') = a(t', x t ', h) 

, In Z N (t',x t ,,h; J) 
( N (t) = , 

where x t > is defined in Eq. (5). The characteristic function of fx can be written 

as 

ft 



r - d 

Ee iuf N (t,x,h) = E e iuf N (o,x ,h) +iuE / e iuf N (t',x t ,,h) f N (t',x t >,h)dt'.{Al) 

Jo dt 



Since 



^0 = i(i-*>' 

d . 1 
dt 



j ^ ^ ^ 



one finds through integration by parts that 



E { e<MCw(t,) ^Ctf(0} = ^ {^^[1 - n*(<&) - 2g (1 - ^(912))]} (42) 

£ {e l <-(*')[^( gi 2 2 ) - 2gJV(9i 2 ) - N' 1 ]} . 



i u 

4 



AN 
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By using (42) in Eq. (41), one finds 

,2^2 ft 



v 2 ri 2 r 

Jj] e iuf N (t,x,h) _ e iuf N (0,x o ,h) \ ^ £ / e »«/W(*') dt' U£\ 

1 Jo 

Ee iuf N (t>) (Q t ,(^ 2 2 )-1) dt' 

Ee iufN ^Q tl (^ 2 )dt'. 



AN 

% u 

Jo 

At t = 0, all sites are decoupled and the central limit theorem for i.i.d. random 
variables implies that 

MO, x, h) jV(0, <x 2 (0, h)). (44) 
The last two terms in Eq. (43) clearly vanish for iV — > oo. For instance, 



AT- 1 



since 

-2 



MV(&) = 0(1) 

for i < i c . Therefore, Eq. (43) yields the following linear integral equation for 
the characteristic function: 

v 2 a 2 /"* 

J? c iuf N (t,x,h) _ g c %uf N (0,x ,h) \ H £ I e iuf N (t',x t ,,h) fal _|_ Q /-^\ 

1 Jo 

whose solution is, keeping into account the initial condition (44), 

2 

e iuf N (t,x,h) _ e -\a 2 {t,x,h) _|_ Q (\y 

□ 

Before concluding this Section, we wish to note that from Eq. (43) one can 
also obtain in a very simple way a well known result for free energy fluctuations 
at zero external field and t < 1 [4, 6], i.e., 

V % = In Z N (t)-N ^2 + ^ Y t - Un (45) 

where Y t is a centered Gaussian random variable of variance 

Indeed, setting -u = \/iVs and x = h = in Eq. (43), one obtains the equation 
Ee ^v t N = l _^. f Ee isr >N (VL t ,{g 2 ) -1) dt' - ^ f Ee is <Q t ,(g 2 )dt'. (46) 
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Since Theorem 1 implies, for vanishing external field and t < 1, 



E (m 2 n) - (en)Y = (£fU 3 2 4> - <£i 2 2 > 2 = 0(1), 




1 -t' 



1 



dt' + o(l), 



from which the result (45) easily follows. 

8 Conclusions and outlook 

We have employed the cavity method to prove a central limit theorem for the 
fluctuations of overlaps and free energy, in a region above the Almeida-Thouless 
line. The key ingredient was provided by the control of the coupled two replica 
system. The open question remains to understand whether and how our method 
can be extended to the entire physically expected high temperature region. 

In the case of vanishing external field, our method can be employed to obtain 
very detailed informations on the system in proximity of the critical point. In 
particular, one can obtain lower and upper bounds on the overlap fluctuations, 
at (5 — 1. We plan to report soon on this [24]. 
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